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Abstract — This study investigates the effects of variations in the fugacities of oxygen and sulfur on 
the partitioning of first series transition metals (V, Cr, Mn, Fe, Co, Ni, and Cu) and W among coexisting 
sulfide melt, silicate melt, and olivine. Experiments were performed at I atm pressure, 1350°C, with the 
fugacities of oxygen and sulfur controlled by mixing C0 2 , CO, and S0 2 gases. Starting compositions 
consisted of a Ca0-Mg0-Al 2 0 3 -Si0 2 -Fe0-Na 2 0 analog for a barred olivine chondrule from an ordinary 
chondrite and a synthetic komatiite. The foJh 2 conditions ranged from log /<> 2 = -7.9 to -10.6, with 
log / s values ranging from - 1.0 to -2.5. Our experimental results demonstrate that the f 0 J fs 2 dependen- 
cies of sulfide melt/silicate melt partition coefficients for the first series transition metals are proportional 
to their valence states. The fojfs , dependencies for the partitioning of Fe, Co, Ni, and Cu are weaker 
than predicted on the basis of their valence states. Variations in foj fs 2 conditions have no significant 
effect on olivine/melt partitioning other than those resulting from /o 2 -induced changes in the valence 
state of a given element. The strong foJfs 2 dependence for the olivine/silicate melt partitioning of V is 
attributable to a change of valence state, from 4+ to 3+, with decreasing f 0l . Our experimentally 
determined partition coefficients arc used to develop models for the segregation of sulfide and metal 
from the silicate portion of the early Earth and the Shergottite parent body (Mars). We find that the 
influence of S is not sufficient to explain the overabundance of siderophile and chalcophile elements 
that remained in the mantle of the Earth following core formation. Important constraints on core formation 
in Mars are provided by. our experimental determination of the partitioning of Cu between silicate and 
sulfide melts. When combined with existing estimates for siderophile element abundances in the Martian 
mantle and a mass balance constraint from Fe, the experiments allow a determination of the mass of 
the Martian core ( — 17 to 22 wt% of the planet) and its S content (—0.4 wt%). These modeling results 
indicate that Mars is depleted in S, and that its core is solid. Copyright © 1997 Elsevier Science Ltd 


I. INTRODUCTION 

The segregation of metallic cores from silicate mantles is 
the most significant differentiation event in the evolution of 
the terrestrial planets. Tn the Earth, for example, the inner 
and outer cores combined comprise —33% of the mass of 
the planet. In contrast, the extraction of basalt from the upper 
mantle at mid-ocean ridges produces an oceanic crust com- 
prising only —0.08% of the mass of the Earth. The absolute 
and relative abundances of siderophile elements (elements 
that partition strongly into metallic phases relative to sili- 
cates) that remain in the silicate mantle of a planet following 
core formation provide the only direct geochemical evidence 
relating to the segregation process. These abundances 
should, in principle, reflect equilibration between metallic 
phases and mantle silicates during segregation of material 
into the core. 

The conditions under which core formation took place 
could have had a significant effect on the present-day sidero- 
phile element abundance pattern of the mantle through varia- 
tions in metal/silicate partitioning. The effects of such vari- 
ables as temperature and pressure on siderophile element 
partitioning are beginning to be understood and quantified 
(Jones and Walker, 1991; Murthy, 1991; Keppler and Rubie, 
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1993; Walker et al., 1993; Capobianco and Amelin, 1994; 
Hillgren et al., 1994; Thibault and Walter, 1995; Li and 
Agee, 1996). In this study, we present experiments that 
investigate the effect of variations in the fugacities of oxy- 
gen (/o 2 ) and sulfur (/s 2 ) on the partitioning of moderately 
siderophile elements (V, Cr, Mn, Fe, Co, Ni, Cu, and W) 
among coexisting olivine, silicate melt, and sulfide melt. The 
experimental results are used to constrain the role of sulfide 
in producing the siderophile element abundance pattern in 
the upper mantle of the Earth and the conditions under which 
core formation took place in the parent body of Shergottite 
meteorites (Mars). 

2. EXPERIMENTAL AND ANALYTICAL METHODS 

Experiments were performed at 1 atm pressure in either a Pt- 
wound or Deltech vertical quenching furnace fitted with a 19 mm 
outer diameter fused-silica muffle tube. Fused silica was chosen to 
minimize interaction with the furnace gases, but had the disadvantage 
of needing to be replaced frequently. End-pieces were fabricated 
from fired pyrophyllite and attached to the muffle tube using General 
Electric Red RTV 106 high temperature silicone rubber adhesive 
sealant. The values of fo 2 and / Sj were controlled by mixing C0 2 , 
CO, and S0 2 gases that were fed into the furnace through a tube in 
the upper pyrophyllite plug (Fig. 1). The gases exited through a 
tube in the lower plug and were fed into a Pyrex flask through a 
silica-glass tube wrapped with heating tape, and then vented into a 
fume hood. Temperature in the furnace hotspot was continuously 
monitored using a Pt-Pt y0 Rh i0 thermocouple calibrated against the 
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Fig. 1. Schematic of the furnace configuration used to perform 
the siderophile element partitioning experiments reported in this 
study. 


melting points of NaCl, Au, and Pd on the 1968 International Practi- 
cal Temperature Scale (Biggar, 1972). The thermocouple was posi- 
tioned on the outside of the muffle tube to avoid contamination by 
S-rich gases, and the temperature difference between the position of 
the thermocouple and the position of the sample was calibrated using 
a second thermocouple. 

The mixture of C0 2 , CO, and S0 2 gases required for each experi- 
ment was calculated by assuming that the gas behaves as an ideal 
solution at high temperature and low pressure and is given in Table 
I. Expressions relating Gibbs free energy change to equilibrium 
constants for independent reactions involving nine gas species (0 2 , 
CO, C0 2 , COS, CS 2 , S 2 , SO, S0 2 , S0 3 ) were written using thermo- 
chemical data from the JANAF tables (Chase et al., 1985). Mass 
balance equations relating each element (C, O, and S) to the stable 
gas species were then used to determine the gas mixture necessary 
to achieve the desired combination of f 0i and / s ,. This calculation 
procedure gives results that are in agreement with those determined 
using the method of steepest descent, a calculation procedure based 
on free energy minimization that has been used in other experimental 
studies involving C0 2 -C0-S0 2 gas mixing (White et al., 1958; 
Haughton et al., 1974). 

Starting materials consisted of either a Ca0-Mg0-A3 2 0 3 -Si0 2 - 
FeO-Na 2 0 analog for a barred olivine chondrule from an ordinary' 
chondrite (Fo86 of Ehlers et al, 1992) or a synthetic Munro Town- 
ship komatiite (KOM; Composition B of Kinzler and Grove, 1985) 
(Table 2) to which synthetic FeS and various combinations of John- 
son-Matthey high-purity NiO, CuO, Cr 2 0 3 , Co 3 0 4 , V 2 0 5 , or W0 3 
were added (Table 1). FeS was synthesized by melting a 2 gram 
stoichiometric mixture of Fe metal sponge and S chips in an evacu- 
ated silica-glass tube. The tube containing the FeS mixture was 
placed in a Deltech vertical quenching furnace at 960°C; the tempera- 
ture was raised gradually to 1 250°C over a period of 2 h, then quickly 
to 1270°C and held at that temperature for 2 h before quenching the 
charge into water. 

Experiments were performed by packing the sample powder into 
a crucible fabricated from a single crystal of gem quality San Carlos 
olivine (Fo 8s) . 91 ) following the methods described by Ehlers et al. 
(1992). Olivine crucibles are preferable to containers made from 
high-purity alumina because they prevent reaction between the cruci- 
ble and charge during olivine-saturated melting experiments. The 
crucible was fastened to the top of a silica-glass rod using 0.2 mm 
diameter Pt wire and placed into the furnace through the lower end 
of the muffle tube w ith the hotspot temperature set to 1078°C. The 
lower pyrophyllite plug was then used to seal the furnace tube and 
the gases w T ere set to the desired flow' rates. After the gases had been 
flowing for 30 min, the sample was raised into the hotspot. A second 
glass rod w'as placed between the benchtop and the rod holding the 
sample, and the hotspot temperature w-as gradually raised to 1350°C 


over a period of ~60 min. At the end of each experiment the lower 
glass rod was removed and the upper rod slid down, quenching 
the experiment against the lower pyrophyllite plug. This method of 
quenching resulted in little or no quench growth in the silicate glass, 
but the sulfide melt generally consisted of a coarse intergrowth of 
quench phases, a common problem in experiments containing sulfide 
melts (e.g., Jones and Drake, 1983; Jones and Walker, 1991; Peach 
and Mathez, 1993). Evidence for the loss of Cu from the charges 
due to volatilization was found in one experiment performed at 
relatively oxidizing conditions (log f 0l = —7.9; log / Sj = —2.0; gas 
flow rate of 1.29 mL/s). This problem was minimized in other 
experiments through the use of low gas flow rates and deep crucibles. 

Experimental products were analyzed using either a 4- or 5-spec- 
trometer JEOL 733 Superprobe at the Massachusetts Institute of 
Technology. Analytical conditions for determination of the major 
elements in all phases were an accelerating voltage of 15 kV and a 
beam current of 10 nA with maximum peak counting times of 40 
s. A 1 /zm spot size was used for analyses of olivine and sulfide, 
while a 10 pm spot size was used to perform glass analyses. Each 
reported silicate analysis represents a mean of 10 spot analyses. On- 
line data reduction was accomplished using the phi-rho-z correction 
scheme. Trace elements in the silicate phases were analyzed sepa- 
rately from the major elements using a beam current of 200 nA 
and peak counting times ranging from 250-700 s, depending on 
concentration. Elements present in the sulfide melt in trace amounts 
were analyzed using beam currents of 50-200 nA and peak counting 
times of 250-500 s, depending on concentration. The detection limit 
for trace elements was taken to be 3 standard deviations above 
background counts. Analyzing experiments containing V required 
corrections for peak overlap interferences from Ti K/3 on V K „ and V K/? 
on Cr KQ . This was accomplished by determining the intensity of the 
interfering peak on the pure oxide of the interfering element (i.e., 
Ti0 2 and V 2 0 5 ) and then scaling the interference for a given analysis 
to the concentration of the interfering element and subtracting it 
from the measured intensity of the element of interest. 

The olivine/silicate melt partition coefficients for Cu were deter- 
mined by secondary ion mass spectrometry using the Cameca IMS 
3f ion microprobe at Woods Hole Oceanographic Institution. Analy- 
ses were performed using a 2 nA beam of O' ions focused to a spot 
size of 15 pm. Positive secondary ions were collected and counted 
by an electron multiplier. Molecular interferences were excluded by 
energy filtering using a ±10-V energy window' and a -90V offset 
(Shimizu and Hart, 1982), and checked by determining the 6:, Cu/ 
A5 Cu ratio for each spot. The relationships between intensity and 
concentration for olivine and glass were determined using experi- 
ment KOM-ls, in which the concentration of Cu in olivine is high 
enough to be determined by electron microprobe. 

Accurate determinations of the compositions of experimentally 
produced sulfide melts were complicated by the presence of coarse 
intergrowths of quench phases. This problem has been dealt with in 
other studies by performing multiple broad-beam analyses (Peach 
and Mathez, 1993), or by rastering the beam over large areas 
(Walker et al, 1993). Our approach to this problem was to analyze 
each quench phase separately, to use image analysis to determine 
their relative proportions, and then to calculate a weighted mean 
composition for the bulk sulfide melt. The uncertainty associated 
with the quench phase proportions was estimated by performing 
replicate determinations and calculating the standard deviation of 
the mean for each quench phase. 

3. EXPERIMENTAL RESULTS 
3.1. Synthesis Experiments 

Experiments were performed on a synthetic Munro Township 
komatiite composition (KOM) (Table 2) to determine the parti- 
tioning behavior of V, Cr, Mn, Fe, Co, Ni, Cu, and W among 
coexisting olivine, sulfide melt, and ultramafic silicate melt. The 
partitioning behavior of Ni was also experimentally determined for 
the CaO-MgO-Al 2 O r SiO 2 -Fe0-Na 2 0 barred olivine chondrule ana- 
log composition (Fo86; Table 2) studied by Ehlers et al (1992). 
The komatiite provides a reasonable analog for the type of high- 
degree partial melt that might be present if the silicate portion of a 
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TABLE 1. Experimental conditions and phase assemblages for partitioning experiments. 


Experiment 

Duration 

(hr) 

1°8 f o 2 

lo 8 f S2 

Gas Flow Rates (mL/s) 

C0 2 : CO : S0 2 Run Products 

Ni-Doped Experiments 





Fo86-8s 5 

48 

-7.9 

-2.5 

1.53 : 0.41 : 0.11 

GI, Oliv, Slf 

Fo86-5s ! 

49 

-7.9 

-2.0 

0.95 : 0.30 : 0.15 

Gl, Oliv, Slf 

Fo86-9s 3 

48 

-8.0 

-1.5 

0.51 : 0.30: 0.18 

Gl, Oliv, Slf, PrEn 

Fo86-15s 3 

49 

-8.6 

-1.5 

0.39 : 0.39 : 0.08 

Gl, Oliv, Slf, PrEn 

Fo86-14s ! 

48 

-9.4 

-1.5 

0.27 : 0.71 : 0.08 

Gl, Oliv, Slf, PrEn 

F086-IOS 1 

48 

-10.3 

-1.5 

0.08 : 1.44 : 0.12 

Gl, Oliv, Slf, PrEn 

Ni-Cr-Cu-Doped Experiments 




KOM- Is 2 

48 

-7.9 

-1.8 

0.46 : 0.17 : 0.10 

GI, Oliv, Sp, Slf 

KOM-6s 2 

48 

-8.0 

-1.5 

0.51 : 0.30 : 0.18 

GI, Oliv, Sp, Slf 

KOM- 1 5s 3 

72 

-8.8 

-1.5 

0.50 : 0.61 : 0.10 

Gl, Oliv, Sp, Slf 

KOM-4S 2 

48 

-9.2 

-1.0 

0.00 : 0.77 : 0.20 

GI, Oliv, Sp 10 , Slf 

KOM-3s 2 

48 

-10.3 

-1.5 

0.08 : 1.44 : 0.12 

GI, Oliv, Sp 10 , Slf 

KOM- 1 Os 4 

72 

-10.3 

-1.5 

0.10 : 1.66 : 0.13 

Gl, Oliv, Sp 10 , Slf 

KOM- 26s 4 * 9 

72 

-7.9 

-2.0 

0.87 : 0.28 : 0.14 



72 

-10.3 

-1.5 

0.10 : 1.66 : 0.13 

Gl, Oliv, Sp 10 , Slf 

KOM-9s 4 

72 

-10.6 

-1.5 

0.00 : 1.28 : 0.10 

Gl, Oliv, Sp 10 , Slf 

Ni-CoW-Doped Experiments 




KOM- 16s 5 

72 

-7.9 

-1.8 

0.46 : 0.17 : 0.10 

Gl, Oliv, Sp, Slf 

KOM-20S 5 

72 

-8.0 

-1.5 

0.51 : 0.30 : 0.18 

Gl, Oliv, Sp, Slf 

KOM- 13s 6 

72 

-9.4 

-1.5 

0.27 : 0.71 ; 0.08 

Gl, Oliv, Sp 10 , Slf, (Tngst?) 

KOM- 11s 6 

72 

-10.3 

-1.5 

0.08 : 1.44 : 0.12 

Gl, Oliv, Sp 10 , Slf, Tngst 

Ni-V-Doped Experiments 





KOM- 19s 7 

72 

-7.9 

-1.8 

0.46 : 0.17 : 0.10 

Gl, Oliv, Sp, Slf 

KOM-23s 7 

74 

-8.0 

-1.5 

0.51 : 0.30 : 0.18 

Gl, Oliv, Sp, Slf 

KOM-22s 7 

90 

-9.1 

-1.5 

0.44 ; 0.78 : 0.10 

Gl, Oliv, Sp 10 , Slf 

KOM- 18s 8 

72 

-10.3 

-1.5 

0.08 : 1.44 : 0.12 

Gl, Oliv, Sp 10 , Slf 


Abbreviations: G1 = glass, Oliv = olivine; Slf = sulfide; PrEn = protoenstatite; Sp = spinel; Tngst — 
tungstate mineral. 

1 Doped with 1.5 wt% NiO; 1.5 wt% FeS. 

2 Doped with 1.5 wt% NiO; 1.5 wt% FeS; 1.0 wt% CuO. 

5 Doped with 1.5 wt% NiO; 1.5 wt% FeS; 1.0 wt% CuO, 1 wt% Cr 2 0 3 . 

4 Doped with 2.5 wt% NiO; 1.5 wt% FeS; 2.0 wt% CuO, 1 wt% Cr 2 0 3 . 

5 Doped with 1.5 wt% NiO; 1.5 wt% FeS; 1.0 wt% Co 3 0 4 , I wt% W0 3 . 

6 Doped with 2.5 wt% NiO; 1.5 wt% FeS; 1.0 wt% Co 3 0 4 , I wt% W0 3 . 

7 Doped with 1.5 wt% NiO; 1.5 wt% FeS; I wt% V 2 0 3 . 

8 Doped with 3.5 wt% NiO; 1.5 wt% FeS; I wt% V 2 05. 

Reversal experiment. 

,0 SpineI occurs within sulfide melt. 


terrestrial planet were largely molten (e.g., Nisbet and Walker, 
1982). The F 086 experiments allow the partitioning behavior of Ni 
in a S-bearing system to be directly compared with the results of 
Ehlers et al. ( 1992), who performed Ni-partitioning experiments in 
a S-free system using the F 086 composition. Experimental condi- 
tions were log f 02 values ranging from —7.9 ( I log unit more reduced 
than the fayalite-magnetite-quartz, FMQ, oxygen buffer) to -10.6 
(the iron-quartz-fayalite, TQF, oxygen buffer) and Tog f$ 2 values 
from - l.Oto -2.5 (Table 1; Fig. 2). At I350°Cboth starting compo- 
sitions produce silicate melts saturated with an immiscible sulfide 
melt and olivine. In addition, the KOM experiments crystallize 
chromian spinel, w'hile protoenstatite occurs in the F 086 experiments 
performed at low foJfs 2 ratio conditions (Table 1 ). The KOM ex- 
periments containing W that were performed at reducing conditions 
also contain a tungstate mineral. 


Olivine grains in the KOM experiments are subhedral to euhedral 
in morphology, tend to be elongate, and are up to ~600 //m in their 
long dimension (Fig. 3). Spinels occur as small (~5-20 j/m), 
equant, euhedral grains. In the high foJfs 2 ratio experiments spinel 
is present in the silicate melt and included in olivine grains, while 
in the low f Q J f &2 ratio experiments it occurs exclusively within the 
sulfide melt (Table l ). The tungstate mineral occurs as small (~5 
prn), spheroidal grains within the sulfide melt. The morphology of 
the sulfide melt changes dramatically with varying /o 2 //s 2 experi- 
mental conditions. In the high /o 2 // Sj ratio experiments, the sulfide 
melt is present as small ( — 10 pm) to moderate size (—250 pm) 
spheres throughout the upper portion of the charge. In the low 
f 0 Jf H ratio experiments the sulfide melt forms large ( — 1,500 yum), 
lens-shaped bodies that commonly occur at the silicate melt /gas 
interface (Fig. 3). 


TABLE 2. Compositions (wt%) of silicate starting materials. 


Si0 2 TiQ 2 A1 2 0 3 Cr 2 0 3 FeO MnO MgO CaO Na^O l^O NiO Total 


F086 47.94 — 7.13 — 8.94 — 31.40 1.39 3.19 — — 99.99 

KOM 50.14 0.33 10.55 0.41 10.98 0.36 17.21 9.49 0.37 0.07 0.09 100.00 
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Fig. 2. Plot of log /o 2 versus log f Sl showing conditions at which 
siderophile element partitioning experiments were performed (filled 
circles). Diagonal solid lines represent constant log flj — log 
/ o 2 2 values. Vertical dashed lines represent the iron-quartz-fayalite 
(IQF), iron-wustite (IW ), and fayalite- magnetite-quartz (FMQ) ox- 
ygen buffers at 1350°C. Horizontal dashed line represents the iron- 
iron sulfide (FeFeS) sulfur buffer at I350°C. 


Olivine grains in the Fo86 experiments are subhedral to euhedral 
in morphology, and tend to be equant. The largest grains are typically 
40-60 /im in diameter. Sulfide occurs as spherical to sub-spherical 
blebs, up to 1 10 /jm in diameter, distributed throughout the silicate 
melt. Pyroxene grains are subhedral to euhedral and tend to be 
elongate. The largest grains are up to —700 pm in length. In the 
experiments where olivine and pyroxene coexist, there is a distinct 
trend of increasing protoenstatite and decreasing olivine abundance 
with decreasing fojfs , ratio. There is also a decrease in the total 
number of crystals w ith decreasing f D J / Sj ratio. 

The composition of the silicate melt varies systematically with 
changing f 0 J fs 2 ratio for both starting compositions. The abun- 
dances of FeO, NiO, CuO, and CoO all decrease with decreasing 
fojfs, ratio, while the abundance of dissolved S increases (Tables 
3-6). The decrease in FeO produces a corresponding increase in 
the concentration of the other major elements. This compositional 
variation leads to low-Ca pyroxene saturation in the Fo86 experi- 
ments performed at low foJJ s, ratio conditions. 

The composition of the sulfide melt also varies as the f Q J / s , ratio 
decreases. The sulfide melt in the high foJfs 2 ratio experiments is 
Ni-rich, while in the low' foJfs 2 experiments it is Fe-rich (Table 7; 
Fig. 4). The concentrations of V, Cr, and Mn in the sulfide melt 
increase with decreasing fojfs , ratio. In contrast with the other first 
series transition metals^ the abundances of Cu and Co do not vary' 
in a systematic way. In the case of Cu, this may be related to 
the need to use higher doping levels to obtain accurate microprobe 
analyses of experiments performed at low foJfs 2 ratio conditions 
(see footnotes in Table 1 for the doping levels used in each experi- 
ment). The sulfur contents of the sulfide melts range from 28.8 to 
38.3 wt %, and increase systematically with decreasing foJfs 2 ratio 
(Table 7). 

The weight ratio sulfide melt/silicate melt partition coefficients 
Moi/snicau- Mdt - j n sulfide melt]/fwt% M in silicate 

melt]) for all of the elements studied (with the possible exception 
of W) are strong functions of the fojfs 2 ratio (Tables 8-9; Fig. 
5). The weight ratio olivine/silicate partition coefficients 
(poimne/siiHate Men _ m j n olivine 1 /[ wt% M in silicate melt]) 
for Ni, Mn, Cr, and Cu do not show a systematic variation w-ith 
changing fojfs, conditions, while those for V and W increase with 
decreasing foJfs i ratio (Table 10; Fig. 6). The Co partition coeffi- 
cients are sim’ilar to those determined in a S-free systems by Ehlers 
et al. (1992) for the Fo86 composition. 


Olivine/silicate melt partition coefficients for Ni are within 2 
standard deviations of the values predicted by the model of Kinzler 
et al. ( 1990) in the experiments performed at conditions more oxidiz- 
ing than the iron-wUstite solid oxygen buffer (IW), as is the case 
in the S-free system (Ehlers et al., 1992). The olivine/silicate melt 
partition coefficient in the Fo86 experiment performed at conditions 
slightly more reducing than IW is 4.0 ± 0.4, significantly lower than 
the value predicted on the basis of the compositions of the coexisting 
olivine and melt (7.1 ± 0.9), but similar to the values reported by 
Ehlers et al. (1992) at comparable /o 2 conditions. 

3.2. Approach to Equilibrium 

The approach to equilibrium represented by our experiments was 
determined through (1) reproducibility of partition coefficients in 
different experiments performed at the same conditions, (2) compar- 
ison of the olivine/melt ^ C/Mg for olivine that crystallized directly 
from the melt with that for the olivine crucible adjacent to the melt, 
and (3) a reversal experiment. 

Multiple experiments w r ere carried out at log f Ql = -7.9, log 
/s 2 = - 1 .8 (KOM-ls; -16s; -19s), at log / 0j = -8.0, log / Sj = - 1.5 
(KOM-6s; -20s; -23s), and at log f 0i = -10.9, log f Si = -1.5 
(KOM-3s; -10s; -11s; -18s) (Table 1). The sulfide melt/silicate 
melt partition coefficients for Ni measured in KOM-ls (540 ± 10) 
and KOM-16s (510 ± 10) overlap at the 2a level, but that measured 
in KOM-19s (410 ± 20) is low. This apparent discrepancy may 
reflect an underestimation of the uncertainty associated with the 
image analysis procedure used to determine the proportions of the 
quench phases. There is good agreement among the sulfide melt/ 
silicate melt partition coefficients for Ni (550 ± 10; 580 ± 30; 550 
± 20) and Mn (0.06 ± 0.03; 0.07 ± 0.01) measured at log f Qi 
- -8.0, log / Sj = -1.5. The same is true in the log /o 2 = 10.3, 

log f$ 2 = — 1.5 experiments for Ni (3600 ± 200; 4400 ± 400; 3800 
± 400). The sulfide melt/silicate melt partition coefficients for Mn 
in experiments KOM-lOs (1.07 ± 0.05), KOM-1 Is (1.09 ± 0.03), 
and KOM-1 8s (0.94 ± 0.04) are similar, but the Mn partition coef- 
ficient determined from KOM-3s ( 1.49 ± 0.04) is high. The sulfide 
melt/silicate melt partition coefficients for Cr agree in KOM-1 0s 
(4.8 ± 0.2) and KOM-18s (4.9 ± 0.2). There is also agreement 
between KOM-3s (7.2 ± 0.3) and KOM-1 Is (6.3 ± 0.5). Although 
the pairs of partition coefficients are similar, they do not overlap at 
the 2 a level. 

Equilibration of the crucible with the melt is an important prereq- 
uisite for equilibration of the entire charge. Table 11 demonstrates 
that there is good agreement between the compositions of the olivine 
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Silicate Melt 
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Fig. 3. Backscattered electron image showing a cross section 
through a siderophile element partitioning experiment performed at 
fo 2 conditions near IW. Co-existing phases are silicate melt, sulfide 
melt, and olivine. Experiment was performed at 1 atm, 1350°C, in 
a crucible fabricated from San Carlos olivine. Scale bar is 1000 gm. 
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TABLE 3. Electron microprobe analyses of silicate phases from experiments performed on Ni-doped F086 starting composition. 


Experiment 

Phase 

Si0 2 

ai 2 o 3 

FeO* 

MgO 

CaO 

NajO 

NiO 

S 

Total 

Fo86-8s 

G1 

60.9(8) 

11.3(2) 

9.1(4) 

13.5(7) 

2-5(1) 

2.4(2) 

0.130(9) 

0.031(3) 

99.86 


Oliv 

40.5(4) 

0.05(2) 

9.7(2) 

49.2(5) 

0.04(1) 

— 

0.97(4) 

— 

100.46 

Fo86-5s 

G1 

61.6(3) 

11.7(1) 

7.9(2) 

13.2(2) 

2.34(6) 

3.2(3) 

0.116(3) 

0.028(1) 

100.08 


Oliv 

40.7(3) 

0.08(2) 

9.5(1) 

49.5(3) 

0.04(2) 

— 

0.84(2) 

— 

100.66 

Fo86-9s 

G1 

60.8(4) 

11.5(1) 

8.6(1) 

14.6(1) 

2.32(3) 

2.5(1) 

0.091(7) 

0.041(4) 

100.45 


Oliv 

40.5(2) 

0.1(2) 

9.7(1) 

49.8(5) 

0.05(2) 

— 

0.59(2) 

— 

100.74 

Fo86-15s 

G1 

59.5(3) 

13.2(2) 

8.2(2) 

14.4(1) 

2.62(6) 

2.76(6) 

0.040(4) 

0.047(5) 

100.77 


Oliv 

40.7(3) 

0.03(1) 

8.9(1) 

50.6(5) 

0.05(1) 

— 

0.28(2) 

— 

100.56 

Fo86-14s 

G1 

59.6(2) 

14.9(1) 

6.8(1) 

13.6(1) 

3.1(1) 

2.61(7) 

0.017(2) 

0.060(1) 

100.69 


Oliv 

41.3(3) 

0.04(2) 

7.9(1) 

5M(4) 

0.05(1) 

— 

0.120 (5) 

— 

100.51 

Fo86- 1 0s 

G1 

59.1(7) 

16.5(2) 

3.2(2) 

16.0(4) 

3.2(1) 

2.3(1) 

0.007(2) 

0.080(2) 

100.39 


Oliv 

41.8(3) 

0.07(1) 

3.68(7) 

54.9(4) 

0.06(1) 

— 

0.028(3) 

— 

100.54 


Notes: Units in parentheses represent I standard deviation of least units cited on the basis of replicate analyses; thus, 60.9(8) 
should be read as 60.9±0.8. Abbreviations as in Table 1. 


that crystallized from the komatiite melt and the olivine crucible 
within a few //m of the silicate melt. This indicates that although 
Fe/Mg exchange was necessary for the crucible to achieve equilib- 
rium with the silicate melt over the range of f 0 J /s 2 conditions inves- 
tigated, a close approach to equilibrium was achieved during each 
experiment. 

A reversal of equilibrium was attempted by performing an experi- 
ment at log f 0l = -7.9, log fs 2 = -1.5 for 72 h, then changing the 
gas flow rates to achieve log / 0j = —10.3, log /s 2 = — 1.5 and 
continuing the experiment for an additional 72 h. The run products 
from this experiment consist of glass, olivine, sulfide, and spinel. 
The abundance and morphology of the phases are identical to those 
found in synthesis experiments, with the sole exception of the pres- 
ence of spinel inclusions in some of the olivine grains. The presence 
of these inclusions is an indication that their host olivine grew at 
the initial conditions of the experiment, during which spinel occurred 
in the silicate melt and as inclusions in olivine grains. Analyses of 
the cores of these grains (Foy^+o^), which are 110 to 145 gm in 
diameter, are in good agreement with the composition of the rims 
and the crucible adjacent to the glass (Table 1 1 ), indicating reequili- 
bration. The sulfide melt/silicate melt partition coefficients for Cr 
(4.8 ± 0.2 versus 4.9 ± 0.2), Mn (1.07 ± 0.05 versus 0.99 ± 0.05), 
and Cu (900 ± 100 versus 670 ± 70) in KOM-IOs and the reversal 
experiment overlap at the 2 a level. The partition coefficient for Ni 
in the reversal experiment is low (2400 ± 300) relative to that in 
KOM-IOs (3600 ± 200), but they overlap at 2.4a. Once again, this 
may reflect an underestimation of the uncertainty associated with 
determining the proportions of the sulfide quench phases using image 
analysis. 


4. DISCUSSION 

4.1. Effect of Variable fo 2 and /s 2 on Siderophile 
Element Partitioning 

Explicit dependence of sulfide melt/silicate melt parti- 
tioning on fo 2 and / Sj is governed by an exchange reaction 
of the form: 


M x O y i,icale Me, ‘ + — S?* 


Z 

2 


= s ! ulfi<jc Melt 



(I) 


where M is a metal cation and X, T, and Z are stoichiometric 
coefficients. An increase in fs 2 at a constant f Qz value (i.e., 
a decrease in the f 0 J f$ 2 ratio) will drive this reaction to the 
product side, increasing the sulfide melt/silicate melt parti- 
tion coefficient, while an increase in fo 2 at a constant 
will favor the reactants, leading to a decrease in the sulfide 
melt/silicate melt partition coefficient (e.g., Peach and Ma- 
thez, 1993). This relationship can be quantified by rewriting 
Eqn. I in terms of a single metal cation: 


MO*/i ca,c Melt 



MS i/x lde Mel1 



( 2 ) 


TABLE 4 . Electron microprobe analyses of silicate phases from experiments performed on KOM starting composition doped with Cr, Ni, and Cu. 


Experiment Phase 

Si0 2 

Ti0 2 

ai 2 o 3 

Cr 2 0 3 

FeO 

MnO 

MgO 

CaO 

Na^ 

k 2 o 

NiO 

CuO 

S 

Total 

KOM- Is 

G1 

51.4(3) 

0.38(2) 

10.9(1) 

0.246(3) 

11.2(2) 

0.314(2) 

15.2(1) 

10.0(1) 

0.23(4) 

0.06(1) 

0.086(2) 

0.073(3) 

0.065(2) 

100.15 


Oliv 

39.7(3) 

— 

0.07(1) 

0.158(3) 

10.89(9) 

0.266(4) 

47.8(3) 

0.23(2) 

— 

— 

0.58(1) 

0.007(2) 

— 

99.69 

KOM-6s 

G1 

51.5(1) 

0.35(2) 

10.53(8) 

0.371(3) 

11-2(1) 

0.296(2) 

17.19(9) 

9.33(9) 

0.12(2) 

0.02(1) 

0.096(2) 

0.026(2) 

0.081(2) 

101.11 


Oliv 

40.3(1) 

— 

0.07(3) 

0.22(1) 

9.7(1) 

0.225(2) 

49.2(2) 

0.23(1) 

— 

— 

0.551(3) 

— 

— 

100.50 

KOM- 15s 

G1 

51.5(2) 

0.34(2) 

10.47(5) 

0.530(3) 

9.4(1) 

0.300(4) 

17.65(6) 

9.32(6) 

0.08(1) 

0.01(0) 

0.024(2) 

0.011(2) 

0.131(1) 

99.77 


Oliv 

40.4(2) 

— 

0.12(2) 

0.324(5) 

8.44(7) 

0.214(3) 

49.8(3) 

0.21(1) 

— 

— 

0.134(5) 

— 

— 

99.64 

KOM-4s 

G1 

53.6(1) 

0.35(2) 

11.0(1) 

0.271(1) 

6.8(2) 

0.322(4) 

18.4(1) 

9.84(5) 

0.11(2) 

0.01(1) 

0.010(1) 

0.011(2) 

0.140(2) 

100.86 


Oliv 

41.3(2) 

— 

0.04(2) 

0.178(5) 

6.1(3) 

0.230(2) 

52.5(4) 

0.21(1) 

— 

— 

0.062(4) 

— 

— 

100.62 

KOM- 3s 

G1 

55.5(2) 

0.36(2) 

11.4(1) 

0.222(1) 

3.72(6) 

0.311(3) 

18.38(6) 

10.3(1) 

0.17(2) 

0.03(1) 

— 

0.0086(5)' 0.154(3) 

100.56 


Oliv 

41.9(4) 

— 

0.03(1) 

0.140(1) 

3.55(6) 

0.220(2) 

54.6(7) 

0.22(1) 

— 

— 

0.017(1) 

— 


100.68 

KOM-IOs 

G1 

55.5(2) 

0.33(1) 

11.22(6) 

0.672(3) 

4.2(1) 

0.289(2) 

18.3(1) 

9.77(8) 

0.05(1) 

0.00 

0.0043(8) 

0.010(2) 

0.140(3) 

100.49 


Oliv 

41.6(2) 

— 

0.08(2) 

0.429(4) 

3.81(7) 

0.200(2) 

54.4(2) 

0.19(1) 

— 

— 

0.030(6) 

— 

— 

100.74 

KOM-26s 

G1 

55.0(1) 

0.34(2) 

11.54(8) 

0.716(2) 

4.09(8) 

0.301(2) 

18.41(9) 

9.9(1) 

0.07(2) 

0.01(0) 

0.007(2) 

0.012(2) 

0.133(5) 

100.53 


Oliv 

41.7(3) 

— 

0.07(2) 

0.458(5) 

3.70(9) 

0.205(2) 

54.4(3) 

0.22(2) 

— 

— 

0.045(7) 

— 

— 

100.80 

KOM-9s 

G1 

56.1(2) 

0.34(2) 

11.23(4) 

0.561(3) 

2.94(5) 

0.263(2) 

19.3(1) 

9.92(4) 

0.08(2) 

0.01(0) 

— 

0.012(2) 

0.151(2) 

100.91 


Oliv 

42.1(2) 

— 

0.07(2) 

0.354(3) 

2.54(5) 

0.179(2) 

55.7(2) 

0.21(1) 

— 

— 

0.017(2) 

— 

— 

101.17 


Notes: Units in parentheses as in Table 2. Abbreviations as in Table 1 . 
'Determined by ion microprobe. 
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TABLE 5. Electron microprobe analyses of silicate phases from experiments on KOM starting composition doped with Co, Ni, and W. 


Experiment Phase Si0 2 TiQ 2 A1 2 0 3 Cr 2 0 3 FeO MnO MgO CaO Na 2 0 K 2 0 NiO CoO W0 3 S Total 


KOM- 16s G1 49.1(2) 0.35(2) 10.59(7) 0.226(5)11.7(2) 0.279(5) 16.0(1) 9.03(4)0.13(2) 0.02(1) 0.121(3) 0.48(2) 1.03(2) 0.064(2) 99.12 

Oliv 39.8(2) — 0.07(0) 0.13(1) 11.2(2) 0.191(9) 46.2(2) 0.22(2) — — 0.67(9) 0.87(8) 0.016(3) — 99.37 

KOM-20s G1 50.4(1) 0.35(2) 10.67(8) 0.254(2)11.1(1) 0.297(3) 16.9(1) 9.23(9)0.05(2) 0.00(1) 0.079(2) 0.435(2) 0.973(8)0.08(1) 100.82 

Oliv 40.5(2) — 0.09(1) 0.151(2) 9.68(9)0.217(4) 49.0(2) 0.22(1) — — 0.45(1) 0.85(1) 0.052(5) — 101.21 

KOM- 13s G1 51.9(3) 0.34(2) 10.6(1) 0.213(3) 8.2(2) 0.286(2) 17.71(8)9.35(6)0.05(1) 0.01(1) 0.018(2) 0.136(3) 0.71(2) 0.136(9) 99.66 

Oliv 41.0(2) — 0.04(2) 0.134(3) 7.6(1) 0.205(4) 50.9(2) 0.20(3) — — 0.097(2) 0.263(5) 0.028(5) — 100.47 

KOM- 11s G1 55.0(1) 0.37(2) 11.1(3) 0.163(1) 3.7(1) 0.285(2) 18.7(1) 9.8(1) 0.04(2) 0.01(1) 0.0034(9) 0.0357(6) 0.669(7)0.141(2) 100.02 

Oliv 41.9(2) — 0.03(2) 0.107(2) 3.39(6)0.200(2) 54.5(2) 0.19(2) — — 0.024(3) 0.070(2) 0.036(6) — 100.45 


Notes: Units in parentheses as in Table 2. Abbreviations as in Table 1. 


The equilibrium constant for this reaction is given by: 

MSy, 


K e u = 


Sulfide rY! IX 

t J o 2 


Silicale f7J 2X 

(lMOy tX J s 2 

Taking logarithms and rearranging gives: 


( 3 ) 


log 


„ Sulfide 

a *H(X 


Silicate 
« MO, 


■ = log /i ,2Y - log + log *■„ (4) 


If the ratio of the activity coefficients for the M species in 
the silicate and sulfide melts remains constant over the range 
of foj k 2 conditions of interest, it can be combined with the 
equilibrium constant and a molar-to-weight-ratio conversion 
factor to form a single constant term. Then, if Y ~ Z, Eqn. 
4 can be simplified to give the expression: 

log \ [log f'g - log/ <7; I + C (5) 

A 


Therefore, if all of the appropriate assumptions are met, a lin- 
ear relationship will exist between log Dm 1 ™ 6 Mel,/S,1,tfa * e Mclt 
and log / sj 2 — log / o 2 with a slope that is equal to the 
anion-to-cation ratio for the oxide or sulfide species of inter- 
est. The strength of the foJh 2 dependence for a given ele- 
ment should be proportional to its valence state. 

Figure 5 is a plot of log Dm '* 1 ** Mcl,/snka,c McTt versus log 
/ ' s [ 2 — log f o 2 showing the partition coefficients determined 
from our experiments on the synthetic komatiite composi- 
tion. The partition coefficient for W has been omitted be- 
cause it could only be determined for a single experiment. 
The logarithms of the partition coefficients for all of the 
elements studied are strong linear functions of log fl! 2 
— log / o 2 2 , as expected. The slopes of the best-fit regression 
lines steepen from 0.38 ± 0.06 for Cu to 1 .37 ± 0.08 for Cr 


(Table 8). The elements Fe, Co, Ni, and Cu are chalcophile 
(^suffidcM/stoeMeH > j 0) at all of the foj fs 2 conditions 
investigated. The elements V, Cr, and Mn go from being 
incompatible in the sulfide melt at high f Q J fs 2 ratios to being 
compatible at low foJfs 2 conditions. Due to the difference 
in their slopes, the partition coefficients forMn and V reach 
a crossover point at log fs 2 2 - log f o 2 = 4.35, where V 
goes from being less compatible than Mn to being more 
compatible (Fig. 5). There is a similar crossover involving 
Fe and Cr at log f$ 2 — log fU 2 = 5.04. 

The relationship between the fo 2 /fs 2 dependence of sulfide 
melt/silicate melt partitioning and anion-to-cation ratio 
means that any change in the valence state of an element 
will complicate the linear relationship given by Eqn. 5. The 
advantage of partitioning among three coexisting phases is 
that any change in the valence state of an element with 
varying fojfs, conditions inferred from sulfide melt/silicate 
melt partitioning should also be reflected in the olivine/ 
silicate melt partitioning data, providing an independent 
check. Unlike sulfide melt/silicate melt partitioning, olivine/ 
silicate melt partitioning of moderately siderophile elements 
is not governed by an exchange reaction involving S and, 
therefore, should be largely insensitive to variations in / S2 . 
A significant and systematic change in olivine/silicate melt 
partitioning with varying foJfs , conditions should, there- 
fore, indicate a change in valence state related to changing 
/o 2 conditions. The olivine/ silicate melt partition coefficients 
for Cr, Mn, Ni, Co, and Cu do not vary systematically with 
changing foJfa 2 conditions. The small amount of variability 
that does exist may be attributable to changes in the composi- 
tions of the coexisting olivine and silicate melt (e.g.. Hart 
and Davis, 1978; Kinkier el al., 1990; Hirschmann and 
Ghiorso, 1994). Of these four elements, Cr and Mn have 


TABLE 6. Electron microprobe analyses of silicate phases from experiments performed on KOM starting composition doped with V and Ni. 


Experiment Phase Si0 2 Ti0 2 A1 2 0 3 

Cr 2 0 3 

FeO 

MnO MgO 

CaO 

Na 2 0 

K 2 0 

NiO 

v 2 o 5 

S 

Total 

KOM- 19s 

G1 

50.0(2) 0.346(3) 10.6(1) 

0.240(2) 

11.4(2) 

0.288(3) 16.6(1) 

9.37(7) 

0.14(3) 

0.03(1) 

0.176(3) 

1.00(1) 

0.062(3) 

100.25 


Oliv 

40.1(3) 0.005(1) 0.07(0) 

0.148(3) 

10.0(2) 

0.212(2) 47.9(3) 

0.24(2) 

— 

— 

1.05(2) 

0.096(3) 

— 

99.82 

KOM-23s 

G! 

49.6(2) 0.347(3) 10.64(9) 

0.250(2) 

11.9(2) 

0.293(3) 17.2(1) 

8.81(7) 

0.15(2) 

0.02(1) 

0.095(3) 

1.048(5) 

0.101(6) 

100.45 


Oliv 

40.4(2) 0.004(1) 0.06(2) 

0.145(2) 

10.2(1) 

0.217(2) 48.7(2) 

0.22(2) 

— 

— 

0.52(2) 

0.104(2) 

— 

100.57 

KOM-22S 

G1 

51.3(2) 0.356(4) 10.7(1) 

0.247(6) 

9.9(2) 

0.303(5) 17.9(1) 

9.17(6) 

0.06(1) 

0.01(1) 0.027(5) 

0.99(1) 

0.145(3) 

101.11 


Oliv 

40.8(2) 0.0047(7) 0.07(2) 

0.147(1) 

8.9(1) 

0.216(3) 50.4(2) 

0.21(1) 

— 

— 

0.152(4) 

0.157(2) 

— 

101.06 

KOM- 18s 

G1 

55.0(2) 0.367(2) 11.15(8) 

0.162(2) 

3.87(7) 

0.290(2) 18.64(8) 

9.9(1) 

0.09(2) 

0.02(1) 0.006(1) 

0.811(5) 

0.153(1) 

100.45 


Oliv 

41.6(3) 0.0062(8) 0.06(2) 

0.105(2) 

3.4(1) 

0.205(2) 54.3(2) 

0.19(2) 

— 

— 

0.027(2) 

0.276(1) 

— 

100.17 


Notes: Units in parentheses as in Table 2. Abbreviations as in Table I. 
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TABLE 7. Electron microprobe analyses of experimentally-produced sulfide melts. 


Experiment 

Fe 

Ni 

Cu 

Cr 

Mn 

Co 

W 

V 

S 

Total 

Ni- Doped Experiments 










F 086 - 8 S 1 

13.3(3) 

58.6(3) 

— 

— 

— 

— 

— 

— 

28.8(2) 

100.7 

Fo86-5s 

14.9(3) 

56.3(8) 

— 

— 

— 

— 

— 

— 

29.6(2) 

100.8 

Fo86-9s 

21.9(6) 

46.5(6) 

— 

— 

— 

— 

— 

— 

31.3(1) 

99.7 

Fo86-15s 

42.2(12) 

24.0(9) 

— 

— 

— 

— 

— 

— 

34.7(3) 

100.9 

Fo86-14s 

43.9(5) 

21.8(4) 

— 

— 

— 

— 

— 

— 

33.9(3) 

99.6 

FO86-10S 1 

50.7(8) 

11 . 2 ( 6 ) 

— 

— 

— 

— 

— 

— 

38.3(5) 

100.2 

Ni-Cr-Cu-Doped Experiments 









KOM-ls 

19.9(2) 

36.4(2) 

14.8(2) 

2 

2 

— 

— 

— 

30.0(1) 

101,10 

KOM- 6 s 

20.9(5) 

41.6(7) 

6.5(1) 

0.034(8) 

0.013(8) 

— 

— 

— 

31.7(3) 

100.75 

KOM-15s 

36.7(7) 

21 . 2 ( 6 ) 

4.8(2) 

0.14(5) 

0.060(6) 

— 

— 

— 

36.6(1) 

99.50 

KOM-4$ 

44.6(6) 

11-2(5) 

6.1(3) 

0.4(1) 

0.17(1) 

— 

— 

— 

36.8(2) 

99.27 

KOM-3s 

51.0(3) 

7.9(2) 

6.5(3) 

1 . 1 ( 1 ) 

0.36(1) 

— 

— 

— 

34.3(2) 

101.16 

KOM-lOs 

43.3(5) 

12.1(3) 

6.8(3) 

2 . 2 ( 2 ) 

0.24(1) 

— 

— 

— 

34.9(2) 

99.54 

KOM-26s 

42.5(7) 

13.2(4) 

6.6(4) 

2.4(3) 

0.23(1) 

— 

— 

— 

34.9(2) 

99.83 

KOM-9s 

43.1(6) 

10 . 0 ( 2 ) 

9.3(3) 

2 . 6 ( 2 ) 

0.37(3) 

— 

— 

— 

34.1(4) 

99.47 

Ni-Co-W-Doped Experiments 









KOM-16s 

13.9(4) 

48.5(4) 

— 

3 

3 

6.7(1) 

2 

— 

30.1(2) 

99.20 

KOM-20s 

23.5(3) 

36.0(5) 

— 

2 

2 

8.3(2) 

2 

— 

33.0(3) 

100.80 

KOM-13 

36.6(4) 

22.1(4) 

— 

0.26(3) 

0.067(9) 

6 . 6 ( 1 ) 

2 

— 

34.7(2) 

100.33 

KOM-lls 

48.6(6) 

11-7(3) 

— 

0.7(2) 

0.24(1) 

4.6(1) 

0.06(1) 

— 

34.7(2) 

100.60 

Ni-V-Doped Experiments 










KOM-19S 

14.0(3) 

56.4(3) 

— 

3 

3 

— 

— 

0.014(1) 

30.4(3) 

100.81 

KOM-23S 

25.9(8) 

41.1(9) 

— 

2 

0.015(6) 

— 

— 

0.013(1) 

33.3(4) 

100.33 

KOM-22s 

39.9(7) 

24.8(6) 

— 

0.18(4) 

0.060(4) 

— 

— 

0.10(2) 

35.9(2) 

100.94 

KOM-18S 

47.5(3) 

17.7(3) 

— 

0.54(6) 

0.21(1) 

— 

— 

0.6(2) 

34.4(2) 

100.95 


Notes: Units in parentheses as in Table 2. Abbreviations as in Table 1 . 
•Sulfide melt quenched to a glass. 

Concentration is at or below detection limit. 

3 Not determined. 


log D s M u,flde MeU/s,1,calc MeU versus log /s' 2 - log / o 2 slopes that 
are consistent with their probable valence states (3+ and 
2+, respectively; Table 8). 

The Cr result is especially significant. The behavior of Cr 
in S-free systems has been investigated in studies of both 
terrestrial and lunar basalt compositions ( Akella et al., 1976; 
Huebner et al., 1976; Schreiber and Haskin, 1976; Barnes, 
1986; Murck and Campbell, 1986; Seifert and Ringwood, 
1988; Roeder and Reynolds, 1991; Mikouchi et aL, 1994). 
One of the fundamental observations to come out of these 
investigations is that the solubility of Cr in spinel-saturated 
silicate melts varies systematically with temperature, melt 


Ni 



Fig. 4. Ternary plot of Fe-Ni-S system comparing compositions 
of sulfide melts from Ni partitioning experiments performed on the 
Fo86 starting composition ( filled circles) with those from the experi- 
mental studies of Rajamani and Naldrett (1978) and Peach and 
Mathez ( 1993) . Shown for reference are the FeS-NiS and FeS-Ni 3 S 2 
joins. 


composition, and fo r Although this effect has been attrib- 
uted primarily to variations in the Cr 2+ /Cr 3+ ratio of the 
silicate melt with changing conditions, the valence state of 
Cr as a function of fo 2 has only been directly measured in 
an Fe-free system (Schreiber and Flaskin, 1976). Our experi- 
mental results are consistent with Cr 3+ being the dominant 
species at 1350°C to conditions as reducing as the IQF oxy- 
gen buffer. 

There are at least three lines of evidence that support a 
low Cr 2+ /Cr 3+ ratio in our silicate melts. The first, men- 
tioned above, is that the slope of log £)gj indcMc1l/s,,,cafc Mcl1 
versus log /s' 2 - log /o' 2 (1.37 ± 0.16) is within 2cr of 
the value expected for Cr 3+ (1.50). Although some percent- 
age (~25%) of the total Cr could be present as Cr 2+ , none 
is required. For comparison, Schreiber and Haskin (1976) 
reported 80% of the total Cr in their silicate melts to be 
divalent in simple system (Ca0-Mg0-Al 2 0 3 -Si0 2 ) experi- 
ments performed at 1500-1550°C and / 0j conditions near 
IQF. The second important observation is that the olivine/ 
silicate melt partition coefficient for Cr measured in our 
experiments (0.62 ± 0.03) does not vary systematically with 
decreasing /> 2 , as might be expected if the Cr 2+ /Cr 3+ ratio 
were changing. A significant increase in the Cr 2+ /Cr 3+ ratio 
of the silicate melt with decreasing />, should result in an 
increase in the olivine/silicate melt partition coefficient 
(Schreiber and Haskin, 1976). Mikouchi et al. (1994) also 
reported a nearly constant olivine /silicate melt partition co- 
efficient for Cr in 1 atm gas-mixing experiments performed 
at 1225 and 1400°C, with /o 2 conditions ranging from 4 log 
units above to 2 log units below the IW buffer. One caveat to 
this observation comes from preliminary experiments which 
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TABLE 8. Regression parameters for logarithm of sulfide melt/silicate melt partition coefficient versus 
one -half the logarithm of the sulfur fugacity/oxygen fugacity ratio. 


Element 

Number of 
Expts. 

. _ Sulfide /Silicate 

log D . = m 

m b 

- log — + b 
2 \ 
r 2 

Calculated 

Valence 

Possible 

Valence 

Weight Ratio Partition Coefficients 





V 

4 

1.36(17) 

-5.9(6) 

0.9742 

2.80* 

3+~4+ 

Cr 

10 

1.37(8) 

-5.3(3) 

0.9733 

2.74+ 

2+-3+ 

Mn 

11 

1.09(6) 

-4.7(2) 

0.9737 

2.18+ 

2+ 

Fe 

15 

0.68(3) 

-1.82(11) 

0.9761 

1.36+ 

2+-3+ 

Co 

4 

0.69(5) 

-0.88(20) 

0.9878 

1.38+ 

2+ 

Ni 

13 

0.65(4) 

0.64(15) 

0.9570 

1.30+ 

2+ 

Cu 

7 

0.38(3) 

1.27(10) 

0.9768 

0.76+ 

1+-2+ 

Molar Ratio Partition Coefficients 





V 

4 

1.36(16) 

-5.6(6) 

0.9728 

2.80+ 

3+-4+ 

Cr 

10 

1.38(8) 

-5.0(3) 

0.9753 

2.76+ 

2+-3+ 

Mn 

11 

1.09(6) 

-4.4(2) 

0.9744 

2.18+ 

2+ 

Fe 

15 

0.68(3) 

-1.50(11) 

0.9753 

1.36+ 

2+-3+ 

Co 

4 

0.69(6) 

-0.55(22) 

0.9859 

1.38+ 

2+ 

Ni 

13 

0.65(4) 

0.97(16) 

0.9547 

1.30+ 

2+ 

Cu 

7 

0.38(2) 

1.59(8) 

0.9846 

0.76+ 

1+-2+ 


Notes: Units in parentheses are lo uncertainties on regression parameters. 


suggest that the olivine/silicate melt partition coefficient for 
Cr 2+ (0.73) may be only moderately larger than that for 
Cr 3+ (0.65) in an Fe-free system (Hanson et al., 1996). The 
third line of evidence for the lack of significant Cr 2+ in our 
experiments is that the abundance of Cr in the silicate melt 
does not vary systematically as a function of /o 2 . The experi- 


ments that were doped with higher concentrations of Cr 
have systematically higher Cr in the silicate melt, but the 
abundance does not increase with decreasing /o 2 . A possible 
explanation for the difference between our results and those 
from simple-system experiments is that the presence of Fe 3+ 
in the silicate melt may shift the Cr 3+ /Cr 2+ transition to 


TABLE 9. Experimentally-determined sulfide melt/silicate melt partition coefficients. 


Experiment 

log/jj 2 - log/o 2 2 D v 

D C r 


D f , 

0C. 

0m 

0 C . 

D w 

Ni- Doped Experiments 








Fo86-8s 

2.70 — 

— 

— 

1.88(2) 

— 

572(7) 

— 

— 

Fo86-5s 

2.95 — 

— 

— 

2.43(2) 

— 

619(5) 

— 

— 

Fo86-9s 

3.25 — 

— 

— 

3.29(3) 

— 

649(15) 

— 

— 

Fo86-15s 

3.55 — 

— 

— 

6.63(7) 

— 

760(30) 

— 

— 

Fo86-14s 

3.95 — 

— 

— 

8.25(6) 

— 

1600(60) 

— 

— 

Fo86-10s 

4.40 — 

— 

— 

20.1(4) 

— 

2160(160) 

— 

— 

Ni-Cr-Cu- Doped Experiments 








ROM- Is 

3.05 — 

— 

■ 

2.29(6) 

— 

540(10) 

250(20) 

— 

KOM-6s 

3.25 — 

0.13(3) 

0.06(3) 

2.40(6) 

— 

550(10) 

313(9) 

— 

ROM- 15s 

3.65 — 

0.39(5) 

0.26(4) 

5.0(6) 

— 

1100(100) 

500(100) 

— 

KOM-ls 

4.10 — 

2.2(2) 

0.68(2) 

8.4(1) 

— 

1430(50) 

690(50) 

— 

KOM-3s 

4.40 — 

7.2(3) 

1.49(4) 

17.6(3) 

— 

— 

960(50) 

— 

ROM-lOs 

4.40 — 

4.8(2) 

1.07(5) 

13.3(3) 

— 

3600(200) 

900(100) 

— 

KOM-26S 1 

4.40 — 

4.9(2) 

0.99(5) 

13.4(5) 

— 

2400(300) 

690(70) 

— 

KOM-9s 

4.55 — 

6.8(3) 

1.8(1) 

18.9(5) 

— 

— 

970(50) 

— 

Ni-Co-W-Doped Experiments 








ROM- 16s 

3.05 — 

— 

— 

1.53(4) 

17.7(5) 

510(11) 

— 

— 

KOM-20S 

3.25 — 

— 

— 

2.72(6) 

24.3(9) 

580(30) 

— 

— 

ROM- 13s 

3.95 — 

1.8(1) 

0.30(2) 

5.7(2) 

62(4) 

1600(100) 

— 

— 

KOM-lls 

4.40 

6.3(5) 

1.09(3) 

16.9(4) 

164(7) 

4400(400) 

— 

0.113(8) 

Ni- V- Doped Experiments 








ROM- 19s 

3.05 0.025(1) 

— 

— 

1.58(11) 

— 

410(20) 

— 

— 

KOM-23s 

3.25 0.022(1) 

— 

0.07(1) 

2.80(7) 

— 

550(20) 

— 

— 

ROM 22s 

3.80 0.18(1) 

11(1) 

0,26(1) 

5.2(2) 

— 

1170(80) 

— 

— 

ROM- 18s 

4.40 1.3(1) 

4.9(2) 

0.94(4) 

15.8(7) 

— 

3800(400) 

— 

— 


Notes: Units in parentheses represent 1 a uncertainties of the least units cited calculated by propagating the standard 
deviation of the mean. 

'Reversal experiment. 
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partitioning of Ni (hexagons), Cu (squares), Co (squares with plus- 
ses), Fe (upside-down triangles), Cr (triangles), Mn (circles), and 
V (squares with crosses) on f Q and /s 2 conditions. Error bars are 
1 cr. Where error bars are not shown, they are smaller than symbols. 


more reducing conditions, as suggested by Schreiber and 
Haskin (1976). It is also possible that the presence of S 2- 
dissolved in the silicate melt affects Cr oxidation-reduction 
equilibria. 

The slopes of the log D M inde Mclt/Slhca,c Mclt versus log f\ {2 
— log / o 2 regression lines for Fe (0.68 ± 0.03 vs. 1 .0), Co 


(0.69 ± 0.05 vs. 1.0), Ni (0.65 ± 0.04 vs. 1.0), and Cu 
(0.38 ± 0.03 vs. 0.5) are shallower than expected in Fig. 5, 
demonstrating that the /o,//$ 2 dependencies of their sulfide 
melt/silicate melt partition coefficients are weaker than 
would be predicted on the basis of their valence states. Peach 
and Mathez (1993) also noted weaker than expected 
fo 2 /fs 2 dependencies for the partitioning of Fe (0.78) and 
Ni (0.60) between sulfide melt and silicate melt in experi- 
ments performed at 8 kbar, I450°C. If these slopes resulted 
solely from the relationship given by Eqn. 5, it would require 
that a significant percentage of each element be dissolved in 
the silicate melt as a neutral species even at f 0l conditions 
1 log unit below the FMQ buffer. This is highly improbable, 
especially given the constancy of the olivine/silicate melt 
partition coefficients for these elements over the range of 
fo 2 conditions investigated. A second possibility is that the 
ratio of the activity coefficients for these elements in sulfide 
and silicate melts (7M ,,catc /y M lfidc ) is variable due to the vari- 
ation in the major element compositions of the coexist- 
ing liquids. This should lead to curvature in the log 

£ Sulfide Melt/Silicaic Melt ys ] Q g j U2 _ | 0 g j U2 re | a tj 0 nship, hOW- 

ever, and none is evident. 

Another mechanism for producing weaker than expected 
foJfs, dependencies is for some percentage of the Fe, Co, 
Ni, and Cu to be bonded with S 2 ~, rather than O 2- , anions 
in the silicate melt, as suggested for Ni by Peach and Mathez 
( 1993). At conditions more reducing than FMQ, S dissolves 
in silicate melts as S 2 ~ anions by displacing O 2- anions that 
are bonded with metal cations such as Fe 2+ (Fincham and 
Richardson, 1954; Haughton et al., 1974). This solubility 
mechanism should leave some percentage of the strongly 


TABLE 10. Experimentally-determined olivine/silicate melt partition coefficients. 


Experiment 

Iog/J - tag/£ 

Dy 

D C r 

D Mn 

°Fe 

D C o 

% 

D C u 

D w 

Ni-Doped Experiments 









Fo86-8$ 

2.70 

— 

— 

— 

1.07(1) 

— 

7.46(19) 

— 

— 

Fo86-5s 

2.95 

— 

— 

— 

1.20(1) 

— 

7.24(7) 

— 

— 

Fo86-9s 

3.25 

— 

— 

— 

1.13(1) 

— 

6.48(16) 

— 

— 

Fo86-15s 

3.55 

— 

— 

— 

1.09(1) 

— 

7.0(3) 

— 

— 

Fo86-14s 

3.95 

— 

— 

— 

1.16(1) 

— 

7.1(3) 

— 

— 

Fo86-10s 

4.40 

— 

— 

— 

1.14(2) 

— 

4.0(3) 

— 

— 

Ni-Cr-Cu-Doped Experiments 









KOM-ls 

3.05 

— 

0.642(4) 

0.847(4) 

0.972(2) 

— 

6.74(6) 

0.096(9) 

— 

KOM-6s 

3.25 

— 

0.593(8) 

0.760(3) 

0.866(5) 

— 

5.74(4) 

0.075(8) 

— 

KOM-ISs 

3.65 

— 

0.611(3) 

0.713(4) 

0.898(5) 

— 

5.6(1) 

0.133(7) 

— 

KOM-4s 

4.10 

— 

0.657(6) 

0.714(4) 

0.897(16) 

— 

6.2(2) 

0.11(2) 

— 

KOM-3S 

4.40 

— 

0.631(2) 

0.707(3) 

0.954(7) 

— 

— 

0.19(3) 

— 

KOM-lOs 

4.40 

— 

0.64(1) 

0.692(3) 

0.907(8) 

— 

7.0(7) 

0.09(1) 

— 

KOM-26S 1 

4.40 

— 

0.640(2) 

0.681(2) 

0.864(7) 

— 

6.4(6) 

— 

— 

KOM-9s 

4.55 

— 

0.631(2) 

0.681(3) 

0.905(9) 

— 

— 

0.08(1) 

— 

Ni-Co-W-Doped Experiments 









KOM-16s 

3.05 

— 

0.58(1) 

0.69(1) 

0.957(7) 

1.81(6) 

5.5(2) 

— 

0.016(1) 

KOM-20s 

3.25 

— 

0.594(3) 

0.731(5) 

0.872(4) 

1.954(4) 

5.7(1) 

— 

0.053(2) 

KOM-13s 

3.95 

— 

0.629(6) 

0.717(5) 

0.927(7) 

1.93(2) 

5.4(2) 

— 

0.039(2) 

KOM-lls 

4.40 

— 

0.656(4) 

0.702(3) 

0.916(10) 1.96(2) 

7.1(6) 

— 

0.054(3) 

Ni-V-Doped Experiments 









KOM 19s 

3.05 

0.096(1) 

0.617(4) 

0.736(3) 

0.877(7) 

— 

5.97(6) 

— 

— 

KOM-23S 

3.25 

0.099(1) 

0.580(3) 

0.741(3) 

0.857(5) 

— 

5.47(8) 

— 

— 

KOM-22s 

3.80 

0.159(1) 

0.595(5) 

0.713(5) 

0.899(7) 

— 

5.6(3) 

— 

— 

KOM-18s 

4.40 

0.340(1) 

0.648(4) 

0.707(2) 

0.879(10) 

~ 

4.5(3) 

— 

— 


Notes: Units in parentheses represent la uncertainties calculated by propagating the standard deviation of the mean. 
‘Reversal experiment. 
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Fig. 6 . Plot of log £) 0,ivine/Silica,e Mel. versus l 0g f 1/2 _ log f \f 2 

illustrating the constancy of the olivine/silicate melt partition coef- 
ficients for of Ni (hexagons), Co (squares with plusses), Fe (upside- 
down triangles), Mn (circles), Cr (triangles), and Cu (squares) 
with changing fo 2 and fs 2 conditions. The systematic increase in the 
olivine/silicate melt partition coefficient for V (squares with 
crosses), and possibly W, (pentagons) is attributable to changing 
valence state with decreasing f 0y Error bars are 1 a. Where error 
bars are not shown, they are smaller than symbols. 


chalcophile elements to bond with S 2 ~ rather than 0 2 ~ 
anions in the silicate melt. If the sulfide melt/silicate melt 
exchange reaction given in Eqn. 1 is modified to take this 
into account, the governing reaction becomes: 

M x o,sl; lica,eMc " + y s?“ 

= MxSlf idcM<:l ‘ + y 05“ + y s5 as (6) 

where Y + Z, = X for divalent cations. Following steps 
similar to Eqns. 2 through 4, the expression relating 

^SulfideMelt/SiHcaleMel. tQ and becomes: 

log Meli/Silicate Melt 

= [(Z 2 - Z, ) log 4 - Y log / 0 J + C (7) 

Different values of Z { can be substituted into this expres- 
sion to account for various percentages of MS dissolved 
in the silicate melt. Experimentally determined log 
£> Sulfide Melt /silicate Meit va | ues can then regressed against (Z 2 

- Zi) log /s, - Y log fo 2 and the slope compared with a 
theoretical value. If, for example, Ni is dissolved in the 
sulfide melt as Ni 3 S 2 (Fig. 4) and the Ni 2+ dissolved in the 
silicate melt is a mixture of NiO and NiS, then substituting 
the correct value for Z { into Eqn. 7 will result in a log 

£ Sulfide Melt /Silicate Melt y§ (Z 2 -Z, ) log /s 2 - Y log fo 2 VQgrCS- 

sion line with a slope of 0. 1 7. 

This approach was used to estimate that —31% of the 
total Ni in the silicate melt must be bonded to S 2_ anions 
in order to explain our /)^ lfldeMc!t/SiHcateMelt values. A similar 


treatment of the sulfide melt/ silicate melt partition coeffi- 
cients for Ni determined experimentally by Peach and 
Mathez (1993) results in -16% of the total Ni in their 
silicate melts being present as NiS. Calculations for the Cu 
and Co data, assuming that all Cu is monovalent and all Co 
is divalent in both liquids, indicates that —20% of the total 
Cu in the silicate melt must be Cu 2 S, and that —31% of the 
Co must be CoS. 

Although these percentages are consistent with the con- 
centration of S dissolved in the silicate melts, there are two 
observations that argue against this mechanism for producing 
the observed foJ fs 2 dependencies. The first is that the pres- 
ence of sulfide species in the silicate melt at the concentra- 
tions required to explain our data should produce a noticeable 
decrease in olivine/silicate melt partition coefficients, and 
none is observed. The second is that the low concentra- 
tion of S in the silicate melts precludes this mechanism as 
an explanation for the behavior of Fe. That the log 

£ Sulfide Melt/ Silicate Melt y& j Q g jU2 _ | Q g J \/ '2 s ] G p es f Qr pe, 

Co, and Ni are within uncertainty of each other suggests 
that a single mechanism may be responsible for all three. 
Capobianco and Amelin (1994) determined metal/silicate 
melt partition coefficients for Ni and Co in the CaO-MgO- 
AI 2 0 3 -Si0 2 system over a range of temperature and fo 2 con- 
ditions. Their experimental results show an /o 2 dependence 
for partitioning at 1300-1425°C that is most readily ex- 
plained by /o 2 -induced changes in silicate melt structure. It 
is possible that the low foJ fs 2 dependencies for Fe, Co, Ni, 
and Cu found in our experiments are related to similar melt 
structural changes. 

A comparison of the fo 2 ffs 2 dependencies of the sulfide 
melt/silicate melt (Fig. 5) and of the olivine/silicate melt 
(Fig. 6) partitioning of V demonstrates that its behavior is 
more complicated than the other elements studied. The 
strong, linear increase in log D? livine/Silicat£ Melt with decreas- 
ing / Qj //s 2 ratio (Fig. 6) indicates that V is undergoing a 
change in valence state related to changing /o 2 , most likely 
from 4+ to 3 4-. This should produce a change in slope in 
the V partitioning data shown in Fig. 5 (cf. Holzheid et al., 
1994; Borisov and Palme, 1995), but none is apparent. It is 


TABLE ll. Forsterite contents and Fe/Mg exchange K D values for San Carlos 
olivine crucibles and experimentally-produced olivines. 


Experiment 

Crucibles 

Fo 

Kd 

Olivines 

Fo 


KOM-ls 

88.7 

0.31 

88.7 

0.31 

KOM-6s 

89.9 

0.31 

90.0 

0.31 

KOM- 15s 

91.2 

0.32 

91.3 

0.32 

KOM4s 

93.9 

0.31 

93.9 

0.31 

KOM-3s 

96.2 

0.35 

96.5 

0.32 

KOM 10s 

96.0 

0.33 

96.2 

0.31 

KOM26s 

96.2 

0.31 

96.3 

0.31 

KOM-9s 

97.4 

0.31 

97.5 

0.30 

KOM- 16s 

88.2 

0.33 

88.0 

0.33 

KOM-20s 

89.9 

0.30 

90.0 

0.30 

KOM- 13s 

92.3 

0.32 

92.3 

0.32 

KOM-lls 

96.5 

0.33 

96.6 

0.32 

KOM- 19s 

89.2 

0.31 

89.5 

0.31 

KOM-23s 

89.6 

0.30 

89.5 

0.30 

KOM-22s 

91.1 

0.31 

91.0 

0.32 

KOM- 18s 

96.5 

0.31 

96.6 

0.30 
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Fig. 7. Plot Of log /)Suind,Nfol</.SilicateMen yersus j Qg j 1 17 _ j Qg 
/ o 2 2 comparing the linear relationship found in the KOM experi- 
ments (open circles) with the curvature found in the experiments 
performed on the F086 starting composition (squares). The F086 
experiments represented by open squares are olivine saturated, while 
those represented as filled squares are saturated with an olivine 
+ protoenstatite assemblage. Error bars are Irx. Where error bars are 
not shown, they are smaller than symbols. 


probable that additional experiments are needed to better 
define the nature of the relationship between log 

£ Sulfide Melt /Silicate Melt an(J ]()g f 1 /2 _ j og fU2 

The systematics of sulfide melt/silicate melt partitioning 
for Ni determined in the experiments performed on the F086 
composition differ from those found in the KOM experi- 
ments. Figure 7 shows that whereas the Ni data from the 
komatiite experiments define a linear relationship between 

log £)SuIfid, Md./smcatc Mel, apd ]Qg f 1 12 _ | Qg f V2 * therc fe 

curvature in the trend defined by the F086 experiments. The 
sulfide melt /silicate melt partition coefficients from the F086 
experiments are larger than those from the KOM experi- 
ments at high foJfs 2 ratio conditions, and comparable or 
smaller at lower fo 2 lfs 2 ratios. The simplest explanation for 
this is a variation in the ratio of the activity coefficients for 
Ni in silicate and sulfide melts ( y n? , calc / y Nj iride ) with chang- 
ing f 0 J fs 2 conditions. Given that temperature and pressure 
are constant in these experiments, the variation in the activity 
coefficient ratio must be attributable to the compositional 
changes occurring in one or both of the coexisting liquids. 

The high concentrations of Ni in the sulfide melts ( 1 1.2 to 
58.6 wt%) could lead to a non-Ffenrian activity-composition 
relationship, but the V-doped komatiite experiments have a 
similar range in Ni ( 1 7.7 wt% to 56.4 wt%) and do not show 
the same partitioning behavior. Therefore, it is probable that 
the curvature in the Fo86 data is related to the compositional 
changes occurring in the silicate melt. The most significant 
variation in the silicate melt composition is a decrease in the 
concentration of FeO (from 9.1 to 3.68 wt%) with decreas- 
ing /o 2 //s, ratio. Although an even larger variation occurs 
in the FeO contents of the silicate melts from the KOM 
experiments ( 1 1.9 to 2.94 wt% FeO), there is an important 
difference between the two sets of experiments. The Fo86 
experiments performed at high f D J f Si ratio conditions are 


saturated with olivine, while the low foJfs 2 ratio experi- 
ments are saturated with coexisting olivine and low-Ca py- 
roxene. Tli is demonstrates that as the FeO content of the 
silicate melt decreases the activity of Si0 2 increases until 
pyroxene saturation is reached. 

Variations in the sulfide capacity (C s ) of the silicate melt 
(C s = [wt% S dissolved in the silicate meltl X [/o 2 2 / 
f 1?]) can provide information on whether the compositional 
changes in the Fo86 silicate melts are significant enough to 
produce the curvature in Fig. 7. Sulfide capacity is a strong 
function of the activity of the metal cations that bond with 
S 2 in the silicate melt, especially Fe 2+ (Fincham and Rich- 
ardson, 1954; Haughton et a!., 1974), and should indicate 
how the activity of Ni varies in the silicate melt. In binary 
systems C s also decreases with increasing Si0 2 activity, as 
O 2 anions bonded solely to metal cations are replaced by 
O 2- anions in tetrahedral coordination with Si 4 * (Fincham 
and Richardson, 1954; Abraham et al., 1960; Abraham and 
Richardson, 1960). Figure 8 shows log C s vs. FeO content 
of the silicate melt in our experiments, a relationship that 
should be linear. The curvature in the Fo86 data (squares) 
that is apparent in Fig. 8 is independent of any compositional 
changes occurring in the sulfide melt and is, therefore, indic- 
ative of significant changes in the activities of transition 
metal cations and/or Si0 2 in the silicate melt with decreas- 
ing foJf H ratio. 

Although a decrease in C, with increasing Si0 2 activity 
is consistent with other experimental studies, the decrease 
in Moit /silicate Mdi * s contrar y lo w b a t would be expected 

on the basis of the partitioning of Ni between olivine and 
silicate melt (e.g., Hart and Davis, 1978; Kinzler et al., 
1990). There is an inverse relationship between MgO con- 
tent of the melt and D ^ vinc,saiKMC McIt but silicate melts with 
different Si0 2 and the same MgO content do not show a 
significant variations in DS! ivine/Silitfate Me,t - The filled squares 



Fig. 8. Plot of the logarithm of the sulfide capacity of the silicate 
melt (log C s ) versus FeO content comparing the linear relationship 
found in the KOM experiments (open circles) with the curvature 
found in the experiments performed on the Fo86 starting composition 
(squares). The Fo86 experiments represented by open squares are 
olivine saturated, while those represented as filled squares are satu- 
rated with an olivine -l- protoenstatite assemblage, ler error bars are 
smaller than symbols. 
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in Figs. 7 and 8 represent Fo86 experiments in which the 
silicate melt is saturated with coexisting olivine and protoen- 
statite. The activity of Si0 2 in these experiments is buffered 
by the coexisting crystalline phases. The relationship be- 
tween log C, and FeO for this subset of experiments is non- 
linear, indicating that the activity of transition metals in the 
silicate melt varies substantially even though the activity of 
Si0 2 does not. Therefore, the variation in Ni activity must be 
related to some other compositional variable. One important 
difference between the Fo86 and KOM silicate melts is that 
the former have unusually low CaO relative to naturally 
occurring basaltic or ultramafic silicate melts, causing some 
of the MgO to form MgAl 2 0 4 complexes in order to stabilize 
A1 in the aluminosilicate network (e.g., Bottinga and Weill, 
1972). A spectroscopic study of aluminosilicate glasses 
along the Si0 2 -MgAl 2 0 4 join found A1 in fivefold and sixfold 
coordination with oxygen rather than the tetrahedral coordi- 
nation that characterizes alkali aluminosilicate glasses 
(McMillan and Kirkpatrick, 1992), suggesting that the pres- 
ence of MgAI 2 0 4 complexes may influence silicate melt 
structure. Good linear correlations exist between the calcu- 
lated mole fraction of MgAl 2 0 4 in the silicate melt and log 
£ Sulfide Mei i /silicate Men ^2 _ q. 9047), and between the molar 

Fe0/MgA! 2 0 4 ratio of the silicate melt and log C s (r 2 
= 0.9123) for the F086 experiments. Although correlation 
does not demonstrate cause, it is suggestive that the presence 
of a MgAl 2 0 4 component in the silicate melt may influence 
the activity of Ni. These results demonstrate that, in certain 
cases, experiments performed in simplified analog systems 
may predict more complicated partitioning behavior than is 
found in natural systems. 

4.2, Modeling Core Formation 

In the following sections, core formation in the Earth and 
in the parent body of Shergottite meteorites (Mars) is mod- 
eled as a single-stage equilibrium segregation of metallic 
phases from mantle silicates using the approach developed 
by Rammensee and Wanke ( 1977) as modified by Newsom 
( 1985). Our experimentally determined sulfide melt/silicate 
melt and olivine/silicate melt partition coefficients are com- 
bined with the experimentally determined metal /silicate melt 
partition coefficients of Schmidt et al. ( 1989) and of Drake 
et al. (1989) to evaluate whether the siderophile element 
abundance pattern in the mantle of either planet is consistent 
with the equilibrium segregation of sulfide melt ± metal 
from an olivine-rich mantle at low pressures and tempera- 
tures. To account for the differences in experimental temper- 
atures among this study ( 1350°C) and that of Schmidt et al. 
(1989; 1300°C) and of Drake et al. (1989; 1260°C) partition 
coefficients were calculated relative to the IW oxygen buffer 
using linear equations such as those given in Table 8. No 
other form of temperature correction was applied to the parti- 
tion coefficients. Although phase composition will undoubt- 
edly affect partitioning (e.g., Hillgren et al., 1996), our ex- 
perimental results indicate that for natural mafic and ul- 
tramafic silicate melts compositional effects will be minor 
compared to variations related to /o 2 //s 2 conditions. 

The bulk (sulfide + metal )/total-silicate partition coeffi- 
cient (D m/S ) is given by: 


D nVs = : f (8) 

Fn q + (C 5 ° / C hq )( 1 - F liq ) 

where D t is the (sulfide + metal )/ silicate melt partition 
coefficient for a given siderophile element, F ]iq is the fraction 
of silicate melt divided by the total fraction of silicates, and 
C s °i/C nq is the silicate solid/silicate melt partition coefficient 
for a given siderophile element. This expression is substi- 
tuted into the equation: 


where X is the weight fraction of sulfide segregated from 
the mantle to form the core, and a is the depletion factor 
for a given element. A depletion factor is the ratio of the 
concentrations of a siderophile element (C f ) to a refractory 
Iithophile element (C,) in Cl chondrites divided by the same 
ratio in the silicate portion of a planet (a = [C,7C ; ]ci/[C, / 
Cyl Mamie -t- crusi ) • Equation 9 predicts the mass of the core that 
is required to produce the observed siderophile element/ 
Iithophile element fractionation as a function of the extent 
of silicate partial melting that accompanied core formation. 
Performing this calculation for a series of siderophile ele- 
ments using known depletion factors and experimentally de- 
termined partition coefficients produces a series of curves 
relating core size to extent of silicate partial melting. If the 
abundances of the siderophile elements in the mantle of a 
planet can be explained by a core-forming process for which 
the given partition coefficients are appropriate, all of the 
curves will intersect at a point coinciding with the mass of 
the core of the planet, and with the extent of silicate partial 
melting that accompanied core segregation. 

4.3. The Role of Sulfide in Core Formation in the 
Earth 

The origin of the siderophile element abundance pattern 
in the upper mantle of the Earth is a long-standing problem 
in geochemistry. Ringwood (1966) first noted that the abso- 
lute and relative abundances of the siderophile elements can- 
not be explained by the equilibrium segregation of metal 
from silicate on the basis of known metal /silicate partition 
coefficients. Four models have subsequently been proposed 
to explain the pattern: ( I ) heterogeneous accretion, in which 
a veneer of chondritic material is accreted follow ing segrega- 
tion of the core and mixed into the upper mantle (Morgan 
et al., 1980; Wanke, 1981); (2) inefficient core formation, 
in which some percentage of metal remains in the mantle 
following segregation of the core, and is subsequently oxi- 
dized (Jones and Drake, 1986); (3) equilibrium segregation 
of a sulfide melt from the upper mantle (Brett, 1984); and 
(4) segregation of core-forming material from a magma 
ocean under extreme pressure or temperature conditions 
(Murthy, 1991; Walker et al., 1993; Li and Agee, 1996). 

The equilibrium sulfide segregation model was conceived 
on the basis of the observations that the density of the outer 
core is 10% lower than that of a pure FeNi alloy at compara- 
ble conditions, requiring the presence of a light alloying 
element or elements (Birch, 1961, 1964), and that sulfide 
melt/silicate melt partition coefficients for siderophile ele- 
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TABLE 12. Parameters used in core segregation modeling. 


Depiction Factors 
V Cr Mn 

Fe 

Co 

Ni 

Cu 

Partition Coefficients 
V Cr Mn Fe 

Co Ni 

Cu 

Earth 

2.0 2.3 5.3 



12.8 

14.0 

12.1 

Metal/Silicate Melt 
0.007 0.12 0.004 — 

71 1200 

13 






Sulfide Melt/Silicate Melt 
0.02 0.08 0.03 — 

17 420 

290 






Olivine/Silicate Melt 
0.34 0.61 0.71 — 

1.9 5.9 

0.11 

Shergittite Parent Body 
— 1.0 1.0 2.18 

14.9 

55.7 

46.4 

Metal/Silicate Melt 
— 0.154 0.0055 3.35 

93 1565 

20 






Sulfide Melt/Silicate Melt 
— 4.96 1.09 14.07 140 3171 

864 






Olivine/Silicate Melt 
— 0.61 0.71 0.90 

1.9 5.9 

0.11 






Low-Ca pyroxene/Silicate Melt 
— — — 0.59 — — 







High-Ca pyroxene/Silicate Melt 
_____ o.25 — — 




ments tend to be smaller than metal/silicate melt partition 
coefficients at comparable conditions (Brett, 1984). Brett 
(1984) proposed that core formation took place through the 
segregation of sulfide melt (23 wt% S; 8 wt% 0) from a 
solid mantle. According to this model the siderophile ele- 
ment abundance pattern of the mantle was established by 
equilibrium partitioning between silicate minerals and sul- 
fide melt at low pressures and temperatures, as the melt 
migrated through the shallow mantle via porous flow. At 
greater depth, the melts aggregated to form larger, faster- 
moving bodies that did not maintain chemical equilibrium 
with the surrounding mantle. The S content of the melt de- 
creased through an unspecified reaction in the deep mantle 
before the melt was added to the core. 

Brett (1984) used partition coefficients from several dif- 
ferent low-pressure experimental studies to test the equilib- 
rium sulfide segregation model against upper mantle sidero- 
phile element abundances. These partition coefficients were 
determined from experiments in which the /o,//s 2 conditions 
were not controlled. Here we revisit the model of Brett 
(1984) and test it using f 0 J /^-sensitive partition coeffi- 
cients. The depletion factors used in our modeling were cal- 
culated using the primitive mantle abundances of Jagoutz et 
al. ( 1979) and the Cl abundances of Anders and Grevesse 
(1989), and are normalized to the refractory lithophile ele- 
ment La (Table 12). Brett ( 1984) proposed that upper man- 
tle conditions during core formation were an / S2 equivalent 
to the FeFeS sulfur buffer (log / s , = —4.10 at 1350°C), and 
an fo 2 value equivalent to the TW oxygen buffer (log f Ql 
= -10.20 at 1350°C). These conditions correspond to a log 
fs 2 2 ~ log f'£ value of 3.05. 

The results from our modeling are presented in Fig. 9. 
Because of the incompatible behavior of V, Cr, and Mn with 
respect to sulfide melt at these conditions, their observed 
fractionation from La in the upper mantle of the Earth could 
only be produced by formation of an unreasonably large core 
(91% to 99% of the Earth’s mass). This does not invalidate 
the sulfide model, however, because it has been argued that 


the depletion pattern for these elements is related to their 
relative volatilities and, therefore, was established by accre- 
tionary processes rather than core formation (Drake et al., 
1989). The sulfide model is fairly successful for explaining 
the Co abundance of the mantle requiring a core that com- 
prises 41 wt% of the Earth and an extent of silicate partial 
melting of 100% accompanying core formation. The cross- 
over point for Ni and Cu at these conditions, however, only 
allows for a core that is 3.6% of the mass of the Earth, with 
core formation accompanied by 96% silicate partial melting. 
This core is too small by an order of magnitude. The lack 
of agreement among Ni, Co, and Cu precludes equilibrium 
segregation of a sulfide melt at low temperatures and pres- 
sures as the appropriate model to explain the siderophile 
element abundance pattern of the mantle of the Earth. Choos- 



Fig. 9. Plot of core mass versus extent of silicate partial melting 
that accompanied core segregation showing the results of equilibrium 
sulfide segregation modeling for V, Cr, Mn, Co, Ni, and Cu deple- 
tions in the upper mantle of the Earth. Core formation f Q J f &2 condi- 
tions tested in this model are log / 0 , = -10.2 and log f Sl = —4.1 
at 1 350°C. 
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ing f 0 Jfs, conditions other than those suggested by Brett 
(1984) does not improve the modeling result, because the 
foJA dependence of log 0^ ndeMd,/StoMdl and log 
^sulfide Melt /snicatc Meit are near |y identical (Table 8). Equilibra- 
tion of a core -forming sulfide melt with mantle silicates at 
low pressures and temperatures cannot explain the mantle’s 
nearly chondritic Ni/Co ratio at any reasonable set of 
fo 2 /fs 2 conditions. 

An alternative to segregation of a single sulfide phase to 
form the core is segregation of a mixture of sulfide melt and 
metal, accompanied by partial melting of the silicate mantle 
(Jones and Drake, 1986; Newsom, 1990). This process 
would also have occurred at low temperatures. We tested a 
model of segregation of a mixture consisting of 27 wt% 
sulfide melt (containing 35 wt% S) and 73 wt% metal. This 
combination will produce an outer core with a S content of 
10 wt%, in agreement with current geophysical estimates 
(e.g., Poirier, 1994). The results from this type of model 
are similar to those for extraction of sulfide alone. A large 
core (88% to 100% of the mass of the Earth) is required to 
explain the abundances of V, Cr, and Mn, while a very small 
core (5.2% of the Earth mass) is needed to explain Ni and 
Cu. Once again the abundance of Co in the upper mantle of 
the Earth can be reasonably accounted for (a 29 wt% core), 
assuming no partial melting of the silicate. This result agrees 
with similar calculations performed by Jones and Drake 
( 1986) and by Newsom (1990). Both of these studies con- 
cluded that some percentage of the segregating material must 
be trapped in the mantle in order for this type of core forma- 
tion to be consistent with known siderophile element abun- 
dances. 

Our experimentally determined partition coefficients are 
not appropriate for testing the possibility that sulfide melt 
segregated from a molten upper mantle at extreme pressure- 
temperature conditions to form the core. Experimental stud- 
ies demonstrate that sulfide melt/silicate melt partition coef- 
ficients for some moderately siderophile elements decrease 
with increasing temperature (Walker et al., 1993) and pres- 
sure (Li and Agee, 1996). Further, the experimental results 
of Li and Agee ( 1996 ) suggest a convergence of the partition 
coefficients for Ni and Co with increasing pressure. Although 
these experimental studies indicate that equilibrium sulfide 
melt/silicate melt partitioning may have played a role in 
forming the core of the Earth, neither provides a completely 
satisfactory resolution to the siderophile element overabun- 
dance problem. The experimentally determined partition co- 
efficients of Walker et al. (1993 ) predict mantle Ge/Ni and 
P/Ni ratios that are too high, and the silicate (34.6 wt% 
Si0 2 ; 48.0 wt% MgO) and sulfide (29.7 wt% S) melts pro- 
duced in the experiments of Li and Agee (1996) are poor 
analogs for the mantle (45.2 wt% Si0 2 ; 38.3 wt% MgO) 
and core (< 12 wt% S), of the Earth, respectively. 

4,4. Size and Composition of the Martian Core 

The Shergottites-Nahklites-Chassigny (SNC) meteorites 
are a group of basaltic achondrites whose major and trace 
element compositions (Stolper, 1979; Stolper and McSween, 
1979) and young crystallization ages (McSween et al., 1979; 
Walker et al., 1979) indicate that their parent body was one 


of the terrestrial planets rather than an asteroid. A compari- 
son of Viking measurements of the Martian atmosphere with 
the noble gases trapped in Shergottite impact melt glasses 
provides the strongest evidence that the Shergottite parent 
body (SPB) is Mars (Bogard and Johnson, 1983). The abun- 
dances of siderophile elements in these meteorites, therefore, 
provide clues to the processes involved in the formation of 
the Martian core. 

The size and composition of the Martian core is poorly 
constrained by geophysical measurements due to the large 
uncertainties associated with the planet’s moment of inertia 
(Bills, 1996). Alternative approaches to determining the size 
and composition of the Martian core can be broadly divided 
into cosmochemical and geochemical models. Cosmochemi- 
cal estimates rely on determinations of the mean density of 
Mars (Bills and Ferrari, 1978), theoretical calculations of 
the temperatures at which equilibrium solids condense from 
a gas of solar composition (Grossman, 1972; Lewis, 1972), 
compositional data from the Viking and Mars 5 missions 
(Owen et al., 1977; Surkov, 1977), and the compositions of 
various classes of chondritic meteorites. Geochemical ap- 
proaches are based, at least in part, on chemical analyses of 
SNC meteorites. 

The least complicated of the cosmochemical approaches 
assumes that condensation of equilibrium solids from a gas 
of solar composition took place in a nebula w'ith an adiabatic 
pressure-temperature gradient. Mars is assumed to have ac- 
creted from a restricted annular zone with a lower mean 
temperature than that from w r hich the Earth accreted (Goet- 
tel, 1983). This model predicts a relatively massive core 
(26.4% of the mass of Mars) that is moderately rich in S 
(17.6 wt%) (Table 13). A modified equilibrium condensa- 
tion model in which feeding zones overlap predicts a similar 
core mass (25.7%) and a higher S content (33.8 wt%; 
Weidenschilling, 1976). The high-S core compositions pre- 
dicted by these models result from the assumption that Mars 
accreted from material in a lower temperature portion of the 
nebula and, thus, contains a higher proportion of volatile 
elements and oxidized Fe than the Earth. A second class of 
cosmochemical model uses mixtures of chondritic meteorites 
to produce a bulk composition that is consistent with the 
existing data on Mars (Anderson, 1972; Anders and Owen, 
1977; Morgan and Anders, 1979; Ringwood, 1981). All 
metal and S is then partitioned into the core. These models 
predict a range of core masses (1 1.9 to 19.0 wt%) and S 
contents (3.5 to 18.6 wt%) (Table 13), depending the as- 
sumptions that go into the model and the data used to con- 
strain it. 

Two groups have used the geochemical approach to esti- 
mate the size and composition of the Martian core. They 
derived similar estimates for the siderophile element abun- 
dance pattern of the SPB mantle on the basis of SNC meteor- 
ite compositions. Both groups concluded that the SPB mantle 
is characterized by a Co/Ni ratio that is — 5 — 1 0 times chon- 
dritic, and pronounced depletions in elements with chalco- 
phile affinities. They used these compositional characteris- 
tics to support the existence of a high-S core. Wanke and 
Dreibus (1988) assumed that Mars contains Fe and Ni in 
chondritic proportions relative to Si and estimated a S abun- 
dance on the basis of their heterogeneous accretion model. 
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TABLE 13. Comparison of model Martian core compositions. 



1 

2 

3 

4 

5 

6 

7 

Mantle + crust 
Mg/(Mg+Fe) 

0.79 

0.77 

0.67 

0.77 

0.66 

0.75 

0.72 

Relative Mass 

73.6 

74.3 

88.1 

81.0 

81.8 

78.3 

80.8 

Core 

Fe 

76.2 

60.4 

72.0 

88.1 

63.7 

77.8 

88.1 

Ni 

6.2 

5.8 

9.3 

8.0 

8.2 

8.0 

11.1 

Co 

— 

— 

— 

0.4 

— 

— 

0.3 

Cu 

— 

— 

— 

— 

— 

— 

0.1 

S 

17.6 

33.8 

18.6 

3.5 

9.3 

14.2 

0.4 

O 

— 

— 

— 

— 

18.7 

— 

— 

Relative Mass 

26.4 

25.7 

11.9 

19.0 

18.2 

21.7 

19.2 


1 : Equilibrium condensation from nebula of solar composition assuming adiabatic pressure-temperature 
gradient (Goettel, 1983); 

2: Equilibrium condensation with modified accretion zones (Weidenschilling, 1972); 

3: Mixture of chondritic meteorites constrained by mean density and moment of inertia (Anderson, 1972); 
4: Seven -component chondrite mixture constrained by Viking and Mars 5 isotopic and compositional data 
(Morgan and Anders, 1979); 

5: 30% Orgueil Cl chondrite + 70% high temperature component mixture (Ringwood, 1981); 

6: Fe mass balance constrained by SNC estimate of mantle FeO content (W&nke and Dreibus, 1988); 

7: This study. 


The compositions of SNC meteorites were used to determine 
the FeO content of the Martian mantle, and the Fe in their 
bulk composition was mass balanced between the SPB man- 
tle and core. They partitioned all Ni and S into the core. 
Their model predicts a core mass of 21.7% and a S content 
of 14.2 wt% (Table 13). Treiman et al. (1986) used analyses 
of SNC meteorites to estimate the geochemical characteris- 
tics of the SPB mantle and discussed the broad constraints 
placed on core formation by the relative abundances of Ni 
and Co. Treiman et al. (1987) modeled the SPB mantle 
abundances of Ni, Co, W, P, and Mo to constrain core forma- 
tion. They used a model in which siderophile elements were 
partitioned among four phases (solid metal, liquid metal, 
solid silicate, liquid silicate) without normalizing their con- 
centrations to that of a refractory lithophile element. Treiman 
et al. ( 1987) estimated that the ambient f Ql conditions during 
SPB core formation were I log unit below the TW buffer, 
and predicted a large core (25 to 35 wt%) composed of a 
50:50 mixture of metal and S-rich metallic liquid ( 12.5 wt% 
S in the core). 

An improved estimate for the size and composition of the 
Martian core can be obtained by combining the Fe mass 
balance approach of Wanke and Dreibus ( 1988) with a trace 
element model similar to that used by Treiman et al. ( 1987) . 
A new constraint is provided by our experimental determina- 
tion of the partitioning behavior of Cu, a strongly chalcophile 
element whose abundance in the SPB mantle should be diag- 
nostic of the S content of the core. Depletion factors were 
calculated for the moderately siderophile elements using the 
SPB mantle abundance estimates of Laul et al. ( 1986) (La) 
and of Wanke and Dreibus (1988) (Cr, Mn, Ni, Co, Cu). 
The metal-to-sulfide ratio of the core-forming material and 
the ambient /o 2 //s 2 conditions during core formation were 
taken as independent variables in the model to determine 
whether a reasonable match to the inferred mantle sidero- 
phile element depletions could be achieved. 

First, broad constraints were placed on Martian core for- 
mation by modeling the abundances of Co, Ni, and Cu at the 


most oxidizing conditions compatible with metal segregation 
(the IW buffer). The log / &2 conditions were arbitrarily set 
at -2.0 (at 1350°C) as a starting point. The best match to 
the SPB mantle depletions for Co, Ni, and Cu was achieved 
with a core that is between 20.2 and 25.9%' of the mass of 
Mars, containing only a small amount of S (^0.4 wt%). 
The extent of silicate partial melting accompanying core 
formation inferred from this model is relatively low (0.2 to 
4.9 wt%). Models involving the segregation of S-rich mate- 
rial (14 to 35 wt% S) permit only a small core (5 to 9 wt%) 
and require large extents of silicate partial melting (60 to 
78%) in order to match the Ni and Cu depletions. These 
models do not reproduce the Co depletion. Varying the ambi- 
ent /o 2 //s 2 conditions does not significantly affect this result. 

TTie preliminary estimates for the size and composition 
of the Martian core derived from modeling trace element 
depletions were refined by assuming a chondritic Fe/La ratio 
for Mars and mass-balancing Fe between mantle and core 
using /o, //independent partition coefficients. Mantle/core 
partition coefficients were calculated for each of the elements 
assumed to be in the core (Fe, Ni, Co, and Cu) using Eqn. 
8. The mantle/core partition coefficient for Fe was modified 
by scaling the oli vine/core partition coefficient to account 
for the presence of orthopyroxene and clinopyroxene using 
the SPB mantle mode of Longhi et al. ( 1992) and olivine/ 
pyroxene partition coefficients determined from a Kilboume 
Hole spinel Iherzolite xenolith (Gaetani, 1996). The mantle/ 
core partition coefficients were then used to calculate the 
concentration of each of the elements in the Martian core 
from their mantle abundance. When the core composition 
summed to less than 100%, it was recalculated at a lower 
fo 2 value until the appropriate core composition was obtained 
(Fig. 10; Table 13). The fo 2 conditions during core forma- 
tion indicated by these calculations correspond to 0.53 log 
units below the IW buffer, with a core representing between 
16.7% and 21.6% of the mass of Mars. The extent of silicate 
partial melting accompanying core formation implied by this 
model is 1 .7 to 6.9%. Due to the low S content of the segre- 
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Fig. 10. Plot of core mass versus extent of silicate partial melting 
that accompanied segregation of a low-S core showing the modeling 
results for Co, Ni, and Cu depletions in the Martian mantle. Core 
formation conditions tested in this model are log / 0j 

= - 10.73 and log / s , = -2.0 at 1350°C. 


gating material, the calculations are insensitive to assump- 
tions about fs 2 conditions. 

Although the fojfs , dependence of the partitioning of W 
between sulfide melt and silicate melt could not be deter- 
mined from our experiments, there is enough partitioning 
data available to perform a first-order test of whether the 
depletion of W in the SPR mantle is consistent with our 
proposed core formation model. Using the / 0 ,-de pendent 
metal/silicate melt partition coefficient for W of Schmidt et 
al. (1989) and our determination of the sulfide melt/silicate 
melt partition coefficient near IW, we found that the W 
abundance in the SPB mantle is consistent with the core 
mass inferred from the other elements if the percentage of 
sulfide in the segregating material is increased slightly to 
~5% (1.8% S in the core) and the extent of silicate partial 
melting is increased to ~30 wt%. 

There is no estimate for the V depletion of the SPB mantle 
and the depletion factors for Cr and Mn are both 1, so that 
their abundances cannot be modeled in the same way as Ni, 
Co, Cu, and W. However, the lack of any significant deple- 
tion for these elements in the SPB mantle can be used to test 
our proposed core formation model if Eqn. 9 is rearranged to 
solve for the depletion factor using the mean core size (19.2 
wt%) and melt fraction (4.3%) obtained from our modeling 
results for Fe, Co, Ni, and Cu. The depletion factors for Cr 
and Mn calculated in this way are 1.08 and 1.01, respec- 
tively, consistent with the abundances inferred from SNC 
meteorites. Models involving segregation of sulfide-rich ma- 
terial (50-100% sulfide melt) produce significant depletions 
in both Cr (a =' 1.97-2.88) and Mn (a = 1.18-1.36) in 
the SPB mantle. Therefore, our modeling results demonstrate 
that the segregation of metal -rich. S-poor material at low 
pressures and temperatures to form the Martian core is con- 
sistent with existing estimates for the concentrations of seven 
of the moderately siderophile elements in the SPB mantle. 

Table 13, which is modified after Table 1 of Longhi et 
al. (1992), is a comparison of the SPB core composition 


calculated in this study with previous estimates. Core compo- 
sitions 1 -5 are based on cosmochemical models, while com- 
position 6 is derived from Fe mass-balance. All of the ex- 
isting models predict a Martian core with a significantly 
higher S content than our composition. The low S content 
of the Martian core implied by our modeling results has 
important implications for our understanding of the internal 
structure of Mars and for its early evolution. First, thermal 
models demonstrate that a core consisting dominantly of 
FeNi metal in a planet the size and age of Mars would be 
solid (Schubert and Spohn, 1990). A solid core is consistent 
with the lack of any significant magnetic field associated 
with Mars (Russell, 1978). The low S content of the core 
also implies that Mars is strongly depleted in S. With the 
exception of Morgan and Anders (1979), whose model core 
composition contains only 3.5 wt% S, all of the previous 
studies concluded that Mars is enriched with respect to S 
relative to the Earth. If the solubility of S in a silicate melt 
is taken to be 1500 ppm, the S/La ratio of Mars would be 
only 0.0065 times chondritic. This contrasts sharply with the 
elevated Martian volatile abundances inferred from studies 
of SNC meteorites (Dreibus and Wanke, 1987), but is con- 
sistent with the noble gas depletions observed in the Martian 
atmosphere (Anders and Owen, 1977; Owen et al., 1977). 
There are two possible explanations for a depletion in S. 
The first is that S was efficiently outgassed from the Martian 
interior during the early evolution of the planet. Once out- 
gassed, S may have escaped from Mars along with the noble 
gases during the loss of a primordial atmosphere (Dreibus 
and Wanke, 1985). A second possibility is that the behavior 
of S during condensation of the solar nebula and accretion 
of the terrestrial planets is more complicated than currently 
assumed. 
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